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Abstract 

Magnetization curve and changes of the single-particle excitation spectra by 
magnetic field are calculated for the periodic Anderson model at half-filling 
in infinite spatial dimension by using the exact diagonalization method. It is 
found that the field-induced insulator-to-metal transition occurs at a critical 
field He, which is of the order of the single ion Kondo temperature. The 
transition is of first order, but could be of second order in the infinite system 
size limit. These results are compared with the experiments on the Kondo 
insulator YbBi2. 
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Study of the strongly correlated electron systems is one of the main issues in condensed 
matter physics. The heavy fermion materials are the typical examples of such systems, in 
which the conduction electrons mix with the almost localized 4f or 5f electrons, and form the 
strongly renormalized quasi-particles which have effective mass of 100 to 1000 times larger 
than the bare value. This strong renormalization is mainly due to the local Kondo-type 
processes, which should be suppressed if an energy gap opens at the Fermi level. 

Recently, the materials called Kondo insulators are attracting much interest. SmS 0, 
SmBe i, YbBi2 1, TmSe §, Ce3Bi4Pt3 I, CeNiSn [0, etc. belong to this category. 
These materials are considered to be the intermediate- valence systems [jll, but show the 
activation type increase of the resistivity at low temperatures, leading to insulating ground 
states, instead of the metallic behaviors in the ordinary heavy fermion or valence-fluctuating 
materials. In these Kondo insulators, a gap opens at Fermi level because of the mixing of the 
conduction electrons with the periodic assembly of the localized electrons. The single-site 
Kondo processes, however, are still possible if the gap is comparable to or smaller than the 
Kondo temperature . In addition, the Kondo processes reduce the mixing, so that the gap 
is renormalized to the value of the order of the Kondo temperature. Thus the Kondo effect 
and the renormalization of the gap take place in a self-consistent way. It is observed that 
the specific heat coefficient 7 (extrapolated from the temperatures higher than the gap) is 
enhanced to some extent in the Kondo insulators 0, which indicates an importance of the 
Kondo-type renormalization in these systems. 

To further investigate the characters of the gap, an interesting experiment was done on 
YbB 12 , in which the application of the strong magnetic field of up to 55T destroys the 
gap and the system goes a transition from insulator to metal. The magnetization starts 
to steeply increase at the critical field. These experiments stimulate us to understand how 
the above-mentioned states with strong renormalization will be modified or destroyed by 
the magnetic field from the theoretical point of view. For these purposes, we utilize the 
recently developed ideas and techniques to treat the strongly correlated electron systems 
in the infinite spatial dimensions [l^. We have calculated the single-particle excitation 



spectra to determine the energy gap, which vanishes at certain critical field ifc, indicating 
the insulator-to-metal transition, and the magnetization curves, which exhibits a steep rise 
at He- These results are consistent with the experiments on YbBi2. Furtheremore, we have 
showed that the Van Vleck susceptibility is finite even in insulators, and is enhanced by the 
electron correlation. 

We start with the periodic Anderson model (PAM), 

ija ia 

+ -7^Y.{ftaCia + h.c) + f/^n/i|n/i|, (1) 

where Ef„ = Ef — ah, h = fisH and H denotes the magnetic field. The other notations are 
standard. We neglect the orbital degeneracy and assume the symmetric case Ef = —U/2 
together with half-filling condition to express the simplest model to the Kondo insulators. 
Note that we apply magnetic field only to f-electrons. The reason for this is that the g 
factor for f-electrons, gj, differ from that for conduction electrons, Qc, so that the total 
magnetization is not a conserved quantity. The total susceptibility is given by x = [(fl'/ " 
9c)l9f]Xfi where X/ is the susceptibility in which the magnetic field is applied only to f- 
electrons ||TT|. If we take Qf = Qc and apply magnetic field to both electrons, x vanishes 



for the half-filling at T 0, which apparently contradicts with experiments. We therefore 
calculate only Xf and denote it simply as x the following. In more realistic models, one 
has to take account of the real values of g factors for both electrons, and also the values of 



the magnetic moments in the crystal field split states ||T2|. 

The first term in eq.([^) denotes the kinetic energy of conduction electrons and is scaled 
in such a way that in (i — > oo limit {d denotes the spatial dimensions) the density of states 



(DOS) has the form p(e) = {2/nW)^l - (e/W)^, where W = 2Vdt, if the Bethe lattice 
with the hopping t parameter is assumed. We, however, regard it as a model DOS for the 
real system. Furthermore, this shape is more suitable for the discussion of the insulator- 
metal transition, since it yields the sharp transition in contrast to the Gaussian DOS. We 
also assume in the following that any magnetic order is suppressed because of an insufficient 



nesting or strong frustrations in real systems, so that we study only the paramagnetic phase. 
W = h = ks = 1 is assumed throughout the present paper. 

We replace Ef hj Ef = Ef + U/2 = 0m the Hamiltonian, assuming the electron-hole 
symmetry, and add the term (f//2) J^ia^fio- to compensate the above modification. But the 
latter term is constant because of the electron-hole symmetry even under finite field since 
Hficr = {nf + am)/ 2 and n/ = 1, where m = nji^ — njii defines the magnetization. The 
thermal Green's function for f-electrons is expressed as 

= /*'G«,(l^4-'-E,(«„)' 

where G'j^lk, ien) = {ien — Ef+ah — V^/ («e„ — e^))"^ denotes the unperturbed f-Green's func- 
tion. In d ^ oo, the lattice problem is reduced to solving a generalized Anderson impurity 
embedded in an effective medium self-consistently The unperturbed Green's function 



for this impurity is given by Gfaiitn) ^ = Gfa{i^n) ^ + So-(ie„). We can fit Gfaiitn) ^ by 
that of the impurity Anderson model with finite number of conduction electron levels and 
energy-dependent mixing (£ = 1 ~ A^^): G';j-^*(ze„)~^ = {lEn- E„ + ah-Y.i>^v'j^/{itn-ei„)) 
Accuracy of this fitting is generally very well in spite of the small number of levels. It 



is because not only the mixing V^'s but also the energy levels e^'s are treated as the fitting 
parameters. As a result, bulk quantities can be obtained with good accuracy in the present 
method, although discreteness of the levels refiects in dynamical quantities. E^j above de- 
notes the effective f level, but the fitting always yields E^j = 0. Although their method is 
applicable also to the finite temperatures, we investigate only the zero temperature prop- 
erties. The ground state of this Anderson model was determined by the modified Lanczos 



method and the thermal Green's function Gfa{itn) was calculated in the form of the 
continued fraction by the Lanczos method ||18[. It is known that A^s = 7 is sufficient to 
obtain a good fit to the thermal Green's function by that with finite levels for the Hubbard 
model [|1^] and the same holds true in the present case. A new estimate of the self-energy is 
obtained by So-(ie„) = Gfa-{ien)^^ — G f^{itn)~^ i which completes the self-consistent loop. 
We first present the single-particle excitation spectrum for [/ = 2 and /i = in Fig. 1(a). 



V = 0.5 is used in the following calculations. It should be noted that the bare f-electron 
peaks at ±f//2 are shifted to deeper positions, in accord with the quantum Monte Carlo 
calculation the numerical renormalization group calculation and the self-consistent 
second order perturbation theory |Tl|]. In addition to them, a Kondo resonance peak is 
formed around e = 0, which is split by the mixing with conduction electrons. The size of the 
gap. Eg are plotted in Fig. 2 as a function of U. Eg is about twice the Kondo temperature 
Tk for the single impurity for large U. Here we define Tk by Tk = 1/4x5 where x is the 



magnetic susceptibility obtained by Bethe Ansatz for the impurity Anderson model ||14 . 

It should be noted that the self-energy So-(ie„) is proportional to — ie„ even in the insu- 
lating phase. This is in contrast to the case of the Hubbard model, in which Il„{ien) oc —ten 



in the metallic phase but oc —i/tn in the insulating phase [|T5|. A reason to this difference 
would be that the gap is formed by the mixing in PAM and it exists even for U = Q whereas 
it is purely due to the correlation in the Hubbard model. 

Therefore, the above spectra are well represented by the renormalized mean-field (MF) 
theory (eq.(|^) below with /i = 0), in which in the mixed band dispersion 



(3) 



is replaced with zV'^ , where z is the renormalization factor, = 1 — 9So-(ie„)/(92e„|e„^o- 
For ?7 = 2 we found z = 0.57. The energy gap is given by Eg = E^'^\ek = —W) — E^~\ek = 
+W) = —W + \/W^~^^4zV^ for the symmetric case. It reads Eg ^ 2zV'^/W when V <t:W. 

Next, we apply magnetic field. Fig. 2 displays the magnetization curve of the f- moment 
m = — for U=0, 1, 2 and 3. The dotted curve is the result for the infinite system 
with U = obtained from 

m = |rfe/(e)[p;t(e)-p/i(e)], (4) 

where /(e) is the Fermi function and Pfai,^) denotes the f-component density of states of a 
spin electrons under magnetic field, 

5 



(upper (lower) sign for e > Ef^ (e < Ef„)), 

K = {E,{e) - Ef,)/^{E,{e)-Ef„Y + W\ (6) 

Ecr{e) = e — V^l^e — Efa) and Ef^ = Ef — ah. This formula gives the f component 
susceptibility Xf = 0.76676, which agrees with the initial slope of the curve for [/ = 
in Fig. 3. 

For U = 0, the spectrum consists of only a few sharp peaks due to the smallness of the 
system, so that the magnetization curve shows redundant steps, although its gross feature 
agrees well with eq.^). For finite U, however, the spectra acquire finer structures and mlh) 
becomes smoother, although a redundant two-step jump is still seen at f/ = 1. Beside these 
artifacts, the magnetization curves indicate sharp increases at the critical field h^, which is 
plotted in Fig. 2 as a function of U. The change of the spectra across he are displayed in 
Fig. 1(b) (h'^hc) and (c) (h'^hc), indicating that the up- and the down- spin bands shift 
with each other and the gap closes. Thus the system undergoes a transition from insulator 
to metal at he. 

We also show in Fig. 4 how the gap Eg{h) closes by magnetic field. It is seen that the gap 
closes as Eg{h) = Eg{0) — 2h, but for f/ > (and maybe also for f/ = around h = 0.26), 
Eg jumps to zero before Eg{0) — 2h vanishes, indicating a first order transition. By a naive 
argument we expect that the Kondo renormalization would be suppressed by the magnetic 
field, so that z and the effective mixing increase, hence the gap increases, whereas the up- 
and down-spin bands are relatively shifted by the field, so that the gap reduces. Hence the 
critical field he is determined as a result of the competition of these two effects. The present 
calculation shows that the gap decreases in proportion to —2h, which indicates that z does 
not change but the gap is closed because of the rigid shift of the each spin band. In fact, we 
found that the self-energy scarcely changes its shape for h < h^. 

Whether the transition is of 1st order is a delicate issue. Investigating the details of 
the spectra shown in Figs. 1(b) and (c), we found that a small peak discontinuously moves 
across the chemical potential, e = 0, when h goes through he. Since we have electron-hole 
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symmetry and even number of electrons in total, no state can appear precisely at e = 0. 
Therefore, in the finite system calculation, a continuous change of h across h^. may cause a 
discontinuous shift of peaks across e = 0, yielding a discontinuous closing of a gap. Thus, 
we expect that the transition could become second order in the infinite size limit, but we 
leave it open for further study. 

Sugiyama, et al. ^ analysed their experiment on YbBi2 and concluded that the gap 
decreases slightly faster than linearly with the field (see the insert in Fig. 4). Since there 
are two activation energies in their resistivity data, they introduced a fine structure in the 
model density of states characterized by two gap values 2Ai = 60 K and 2A2 = 180 K. The 
average of the two gives Eg ~ 120K and the critical field He ~ 450 KOe corresponds to 57 
K, using gj = 1.9. Thus they satisfy Eg ~ 2Hc. Our calculations are in good agreement to 
these features at least qualitatively. 

Finally, we show in Fig. 5 the f-component magnetic susceptibility obtained from the 
initial slopes of Fig. 3, together with that multiplied by Eg and z. Note that the susceptibility 
has finite values despite a finite gap opens, since the total magnetization is not conserved 
in the present model. It is sometimes called the Van Vleck susceptibility. Kontani and 
Yamada gave a discussion that it should be enhanced approximately by z~^. Our result 



shows that it is indeed enhanced by the interaction, although the enhancement becomes 
weaker than or E~^ for strong coupling regime, which may be due to the interspin 
exchange interaction. The present result is qualitatively in agreement with the self-consistent 
perturbational result [ pJ]) . 

In summary, we have calculated the magnetization curve and the single-particle exci- 
tation spectra of the periodic Anderson model at half-filling in infinite dimensions in the 
paramagnetic phase, which is a prototype of the Kondo insulators, and found that the 
insulator-to-metal transition takes place by applying magnetic field. The gap decreases as 
Eg{h) = Eg{0) — 2h, and vanishes at the critical field he, which is found to be of the order 
of Tx/2. It indicates that the gap closes for h < in such a way that each spin band is 
shifted rigidly by ±h. The transition is of first order, but could be of second order in the 



infinite system. The above calculation are mostly in accordance with the observations in 
YbBi2, although more realistic model must be used for detailed comparison. Further study 
of the magnetic field effects to the states away from half-filling will also shed light on the 



metamagnetism in metallic heavy fermion systems ||23|| . 
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FIGURES 

FIG. 1. The single-particle spectrum of PAM for ?7 = 2, F = 0.5 and (a) h = 0, (h) h = 0.09 

and (c) h = 0.1. The dotted line in (a) indicates the spectrum of the mean-field theory with 
reduced mixing zV'^ withz = 0.57. 

FIG. 2. The energy gap Eg, the twice of the critical magnetic field 2hc and the impurity Kondo 
temperature Tk by the Bethe Ansatz are plotted against U. 

FIG. 3. The magnetization curve for U = O(circle), l(triangle), 2(square) and 3(diamond). 

FIG. 4. The energy gap for U = O(circle), U = l(triangle), 2(square) and 3(diamond) are 
plotted as a function of the magnetic field. The dotted lines indicate Eg{h = 0) — 2h. The insert 
is the gap 2Ai obtained by Sugiyama, et al. for YbBi2. 

FIG. 5. The f-component magnetic susceptibility (filled circles) and that multiplied by Eg (open 
circles) and by z (square) scaled at f/ = 0. 
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